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Equilibrium long-ranged charge correlations at the surface 
of a conductor coupled to the electromagnetic radiation 
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This study is related to the fluctuation theory of electromagnetic fields, charges and currents. The 
three-dimensional system under consideration is a semi-infinite conductor, modeled by the jellium, 
in vacuum. In previous theoretical studies it was found that the correlation functions of the surface 
charge density on the conductor decay as the inverse cube of the distance at asymptotically large 
distances. The prefactor to this asymptotic decay was obtained in the classical limit and in the 
quantum case without retardation effects. To describe the retarded regime, we study a more general 
problem of the semi-infinite jellium in thermal equilibrium with the radiated electromagnetic field. 
By using Rytov's fluctuational electrodynamics we show that, for both static and time-dependent 
surface charge correlation functions, the inclusion of retardation effects causes the quantum prefactor 
to take its universal static classical form, for any temperature. 

PACS numbers: 05.30.-d, 52.40.Db, 73.20.Mf, 05.40.-a 



I. INTRODUCTION 



Experimental and theoretical investigations of charged 
systems in thermal equilibrium play a key role in the un- 
derstanding of various fundamental properties of solids 
and liquids, in the bulk as well as on the surface. 
There exist two complementary theoretical approaches 
to Coulomb models: one based on the solution of mi- 
croscopic models, the other based on the assumption 
of validity of macroscopic electrodynamics. The micro- 
scopic description is more laborious and complicated but, 
if available, it can reveal restricted applicability of the 
macroscopic theory. On the other hand, the macroscopic 
phenomenology usually allows us to predict, with much 
less effort, basic features of relatively complicated com- 
plex physical systems. 

In this paper, we study the fluctuations of electro- 
magnetic fields, charges and currents in systems formu- 
lated in the three-dimensional (3D) Cartesian space of 
points r — {x,y,z). We shall deal with semi-infinite 
geometries, inhomogeneous say along the first coordi- 
nate X. It will be sometimes useful to denote the re- 
maining two coordinates normal to a; as R = {y,z). 
The physical situation we are interested in is pictured 
in Fig.l. In its most general formulation, the model 
consists of two semi-infinite media with the frequency- 
dependent dielectric functions ei(w) and e2('-^), in the 
half-space Ai = {r = (a;,R);x > 0} and in the comple- 
mentary half-space A2 = {r = (x,R);a; < 0}, respec- 
tively. The interface between the media is localized at 
a; = 0. Although we shall derive basic formulas for this 



eico) 



£l((0) 







X 



t Unite Mixte de Recherche No 8627-CNRS 

*On leave from Institute of Physics, Slovak Academy of Sciences, 

Bratislava; Electronic address: Ladislav.Saniaj@savba.sk 

t Electronic address: [Bernard. Jancovici@th.u-psud.fr| 



FIG. 1: Two semi-infinite media characterized by dielectric 
functions e\{u) and e2(<j-'). 



general case, the analysis of the results will be done for 
the specific case when the half-space Ai is occupied by a 
Coulomb fluid, namely a jellium, and the half-space A2 
is formed by an impenetrable plain hard wall with the 
vacuum permittivity. For a Coulomb fluid composed of 
charged particles, the average particle density is varying 
at microscopic distances (of the order of the correlation 
length) from the interface. Since we shall be interested 
in macroscopic electrodynamic phenomena at distances 
much larger than the microscopic length scale, we can as- 
sume that the Coulomb fluid and the dielectric function 
are homogeneous in the whole half-space Ai. 

The presence of the Coulomb fluid gives rise to a sur- 
face charge density a on the conductor which has to be 
understood as being the microscopic volume charge den- 
sity integrated on some microscopic depth. It is associ- 
ated with the discontinuity of the electric held E at the 
surface a: = of the conductor. At a point r = (0, R) on 
the surface, in Gauss units, 

4^a(R)=ii;+(R)-i?-(R), (1) 

where the upperscript + (— ) means approaching the sur- 



face through the Hmit x ^ 0+ (x — > 0~). The quantity of 
interest is the surface charge density correlation function 
between two points 
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(a(R)a(R')) 
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X ( [E+{R) - E-{ii)] [i?+(R') - i?-(R')] y m 

where (■■•)"'" represents a truncated equihbrium statis- 
tical average, {AB)'^ — {AB) — {A){B), at the inverse 
temperature (3. The system is translationally invariant 
in the R-plane and so the charge correlation function 
depends only on the distance between the points, 

(a(RMR'))^ = ^(|R-R'|). (3) 

It is useful to introduce the Fourier transform 



5'(q)- /"d2i?e"i-^5'(R), 



(4) 



with q = {qyjQz) being the 2D wave vector. 

For classical Coulomb fluids composed of charged par- 
ticles with the instantaneous Coulomb interactions, S{R) 
can be retrieved by simple macroscopic argument based 
on the electrostatic method of images, giving [l| 



/35ci(i?)--^;^, 



(5) 



where the distance R — |R| is much larger than any mi- 
croscopic length scale, like the particle correlation length. 
Since, in the sense of distributions, the 2D Fourier trans- 
form of l/R"^ is — 27rq, the function /35ci(g) has a kink 
singularity at q = 0, 



PScliq) - — g, 



q 



0. 



(6) 



Note the universal form of the asymptotic formulas, inde- 
pendent of the composition of the Coulomb fluid, which 
is related to specific sum rules for systems with Coulomb 
interaction (for a review, see Ref. la). We see that the 
surface charge correlation on the conductor has a long- 
ranged decay. Such behavior is in contrast with a short- 
ranged, usually exponential, decay which occurs in the 
bulk charge density correlations. The same result IS]) has 
been obtained also in the microscopic language Q. The 
macroscopic formulas for the surface charge density cor- 
relation functions on conductors of various shapes were 
derived by Choquard et al [411. The inverse-power law 
behavior of the surface charge density of type (O causes 
the conductor-shape dependence of the dielectric suscep- 
tibility tensor (which relates the average polarization of 
the system to a constant applied electric field, in the lin- 
ear limit) 0,3- Such a phenomenon is predicted by the 
macroscopic laws of electrostatics [yj . 

The extension of the classical result ([S]), or equiv- 
alently ([S]), to a quantum Coulomb fluid was accom- 
plished in Ref. 7. The model under consideration was 



the jellium, i.e. a system of pointlike particles of charge 
e, mass m and bulk number density n, immersed in a 
uniform neutralizing background of charge density —en; 
the background is assumed to have the vacuum dielec- 
tric constant 1. The dynamical properties of the jel- 
lium have a special feature: there is no viscous damping 
of the long-wavelength plasma oscillations for identically 
charged particles. The absence of damping was crucial 
in the derivation of the large-distance behavior of the 
surface charge correlation function using long-wavelength 
collective modes, namely the nondispersive (if the pres- 
sure term is neglected) bulk plasmons with frequency ujp 
and the surface plasmons with frequency lu^, given by 

In the Fourier space, the obtained asymptotic g — > 
result has the non-universal form [7] 
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[2/(a>,) - fJLOp)] 
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q, q-^O, (8a) 

(8b) 



/H^^coth(^ 



According to the correspondence principle, a quan- 
tum system admits the classical treatment in the high- 
temperature region, which corresponds in our case to 
(3h — > 0. In this limit, the function f{uj) = 1 for any 
Lu and the quantum relation (|8aP indeed reduces to the 
classical one ([6]). 

The surface charge correlation function discussed up 
to now was static, i.e. time taken at the two distinct 
points was the same. The generalization to the time- 
dependent correlation function is accomplished by intro- 
ducing a time-dependent surface charge density, defined 
as the Heisenberg operator 

cr(t, R) = exp(ii/t/n)cr(R) cxp(-iift/S) (9) 

with H being the Hamiltonian of the Coulomb system. 
Since the Heisenberg operators at different times do not 
commute, it is useful to introduce the time-dependent 
surface charge correlation function in a symmetrized form 

[a(R)a(R')]t = ^ (a(i, R)a(0, R') + a(0, R')a(t, R))^ . 

(10) 
This symmetrized correlation function still possesses the 
translational invariance in the R-plane, 

[a(R)a(R')]t = S'qu(i,|R-R'|). (11) 

The asymptotic q —^ formula was shown [3, Q to be 



[2f{uJs) cosjujst) - f{u}p) cosjujpt)] 

P^ciu{i,q) '^ -. q- 

q^O 47r 



(12) 



For i = 0, this formula reduces to the static one (|8ap . 
Note that in the classical limit with f{uj) = 1 for any uj, 



the surface charge correlation function is non-universal 
and exhibits a periodic time dependence of type 



(3S,iit,q) 



[2 cos((jJsi) — cos{ujpt)] 



(13) 



The quantum static and time-dependent results ([8a|) 
and ()12|) . respectively, were derived in the nonretarded 
regime where the speed of light c is taken to be infinitely 
large, c = oo. The effects of the retardation take place 
just at large distances we are interested in, and so they 
should be taken into account. To describe the retarded 
regime with the finite value of c, we study a more general 
problem of the semi-infinite Coulomb system in thermal 
equilibrium with the radiated electromagnetic (EM) field. 

We intend to treat both the Coulomb fiuid and the 
radiation like quantum objects. For the static case, a 
substantial simplification arises in the high-temperature 
limit ph -^ 0. First, according to the correspondence 
principle, both matter and radiation can be treated clas- 
sically. Second, the application of the Bohr-van Leeuwen 
theorem [3, ^^ leads to the decoupling between classical 
matter and radiation, and to an effective elimination of 
the magnetic forces in the matter; for a detailed treat- 
ment of this subject, see Ref. [ill This means that the 
matter can be treated as a classical matter, unaffected by 
radiation, where the charges interact only via the instan- 
taneous Coulomb potential. Because of the absence of 
quantum effects for /3h ^ we expect that the classical 
formula ([6]) will be restored in this limit. How the retar- 
dation effects manifest themselves at a finite temperature 
is an open question which is studied in this paper. All 
that has been said in this paragraph does not apply, in 
general, to time-dependent quantities il2!|. In particular, 
in the presence of the radiation, it can happen that the 
classical time-dependent result ([T^ is not reproduced in 
the limit /3S -* 0. 

To deal with the proposed physical problem of such 
complexity, we shall use a macroscopic theory of equi- 
librium thermal fluctuations of EM fiel d, p ublished by 
Rytov [ij] and further developed in Ref. [ij. This theory 
is presented also in the Course of Theoretical Physics by 
Landau and Lifshitz ^15] of which we adopt the notation. 
Although the EM fluctuation theory was extensively ap- 
plied in investigations of thermally excited surface EM 
waves [16J, [l3] , we did not find a study about the present 
topic concerning the long-range decay of surface charge 
correlation functions. 

The main results obtained in this paper are the fol- 
lowing. The quantum static and time-dependent prefac- 



J 



tors of the small-g behaviors ((8a|) and ([T2)) . respectively, 
remain valid at some intermediate distances, where the 
retardation effects do not play any role. In the strict 
large-distance asymptotic limit, for both static and time- 
dependent surface charge correlation functions, the inclu- 
sion of the retardation effects causes the quantum pref- 
actor to take its static classical form, see Eqs. (O or 

The paper is outlined as follows. In Sec. II, we review 
shortly the EM fluctuation theory and derive the formula 
for the surface charge correlation function between two 
semi-infinite media. The analysis of the asymptotic form 
of this formula for the configuration of interest jellium- 
vacuum is the subject of Sec. III. A brief recapitulation 
and concluding remarks are given in Sec. IV. 



II. DERIVATION OF BASIC FORMULAE 

We consider the (3 + l)-dimensional space with spatial 
vectors r and time t. The physical system of interest is 
a medium and an EM field present in it, which are in 
thermodynamic equilibrium. 

The medium is composed of moving charged particles 
which are assumed to be non-relativistic. In the long- 
wavelength scale much larger that the interparticle dis- 
tances in the medium, its isotropic macroscopic charac- 
teristics are the frequency and (possibly) position depen- 
dent dielectric function e{Lu;r) and permeability fi{Lu;r). 
We shall assume that the medium has no magnetic struc- 
ture, i.e. it is not magnetoactive, and /i = 1. 

The matter is coupled to the EM field. The classical 
EM field potentials form a 4- vector (0, A), where 4>{t,r) 
is the scalar potential and A(i, r) is the vector potential. 
In the considered Weyl gauge cj) — 0, the microscopic 
electric and magnetic fields are given by 



E 



19A 



B = curl A. 



(14) 



The elementary excitations of the quantized EM field are 
described by the photon operators Aj {j — x, y, z) which 
are self-conjugate Bose operators. 

Medium and coupled radiation are in thermal equilib- 
rium. The EM fields and inductions are random variables 
which fluctuate around their mean values. These mean 
values are the quantities obeying macroscopic Maxwell's 
equations. The construction of all types of photon 
Green's functions is based on the retarded Green func- 
tion, defined as the tensor 



iDjk{t]v,v') 



{A, {t, r)ife (0, r')-Ak (0, r')i, {t, r)) , t > 0, 



(15) 



0, 



t < 0. 



Here, j,k = x,y,z, Aj{t,r) denotes the vector-potential 
operator in the Heisenberg picture and the angular brack- 
ets represent equilibrium averaging over the Gibbs dis- 
tribution of the whole system. In what follows, we shall 
work in the Fourier space with respect to time. The 
Fourier transform of the retarded Green function reads 

/■OO 

Djk{uj;r,r') ^ dte'^'Djk{t;r,r'). (16) 

JO 

For media with no magnetic structure, the Green func- 
tion tensor possesses the symmetry 



Djk{(^;r,r') = Dkj{uj;r',r) 



(17) 



Within the framework of the fluctuational electrody- 
namics of Rytov [Tj, [ij, \lM J the retarded Green function 
tensor fulfills the differential equation 



E 

1=1 
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dxjdxi 



5u^-5'. 



to 



e(w;r) 



Afc(w;r,r') 



-47rM,fe(5(r - r'). (18) 



Here, in order to simplify the notation, the vector r = 
{x,y,z) is represented as (xi, X2, 2:3). The differential 
equation (fT8|) must be supplemented by certain bound- 
ary conditions. The second space variable r' and the sec- 
ond index k are not involved in the mathematical opera- 
tions on the tensor, so they only act as parameters. The 
boundary conditions are thus formulated with respect to 
the coordinate r and the Green function Dik{uj; r,r') is 
considered as a vector with the components I ~ x,y,z. 
There is an obvious boundary condition of regularity at 
infinity, |r| — j 00. At an interface between two different 
media, the boundary conditions correspond to the macro- 
scopic requirements that the tangential components of 
the fields E and H = B be continuous. Since the electric 
and magnetic fields are related to the vector potential by 
([T^ . the role of the vector components Ei, up to an irrel- 
evant multiplicative constant, is played by the quantity 



i—Dik{LJ;r,r') 



(19) 



and the role of the vector component Hi is played by the 
quantity 



E 



curlyDj7c(w;r,r'). 



(20) 



Here, we use the notation curljj — J2m ^imjd/dxm with 
^imj being the unit antisymmetric pseudo-tensor. In both 
cases (ITO)) and (^0]) . the tangential components, which are 
continuous at the interface, correspond to indices I = y,z. 
The fluctuation-dissipation theorem tells us that the 
fluctuations of random variables can be expressed in 
terms of the corresponding susceptibilities. For the as- 
sumed symmetry (J17p . the theorem implies 



[E,{v)Ek{v% 



UJ 



coth(/3ftw/2) Im Djfc (w; r, r') , 



(21) 



where the spectral distribution of the electric-field fluc- 
tuations [Ej{r)Ek{r')]uj is the Fourier transform in time 
of the symmetrized (truncated) correlation function 

^^E,{t,r)EkiO,r') + Ek{0,v')E,it,r)y . (22) 

Since the studied problem of two semi-infinite media, 
presented in Fig. 1, is translationally invariant in the 
R-plane perpendicular to the x axis, we introduce the 
Fourier transform of the Green function tensor with the 
wave vector q — {qy, qz): 

DA^; r, r') - J ^e-^^<^-^'^D,k{u:, q; x, x'). 

(23) 
At the time being, the dielectric functions ei(w) and 
e2{'^) are taken as general. The Green function tensor 
for simple planar systems was obtained, as the solution of 
the differential equation ()18p supplemented by the men- 
tioned boundary conditions, in a number of papers, see, 
e.g. . Appendix A of Ref. [1^ or, for multilayers, Ref. 
(l8| . We shall not repeat the derivation, but only write 
the final formulae. Let us define for each of the half-space 
regions the inverse length Kj(Lu,q) (j — 1,2) by 

K^j{uj,q) =q^ - -^ej{Lj), ReKj{Lj,q) > 0; (24) 

from two possible solutions for kj we choose the one 
with the positive real part in order to ensure the regu- 
larity of the Green functions at asymptotically large dis- 
tances from the interface (see below) . We shall only need 
the quantity D^xi^^, q', x, x') for which the previously ob- 
tained results can be summarized as follows: 



• If a:,a:' > 0, 



J-^T.T. "~ 



uj^ei 



Six-x')^'-^^ 



Kie2 + K2ei 



.(25) 



• If X < and x' > 0, 

4:'Kh{cq)'^ 



J-^xr, 



K^X — KxX 



LO^ 



K\ti -I- K2ei 



(26) 



• The case a; > and x' < is dcduciblc from Eq. 
(|26|l by using the symmetry relation (fTT]) . 

• If X, x' < 0, considering the 1 <~^- 2 media exchange 
symmetry, we obtain from Eq. (|25p that 



Dx 



-d{x — X ) 



uj^e2 



CJ^e2'«2 

K2ei -I- Kie2 



,(27) 



The symmetrized surface charge correlation function 
pop is expressible in terms of the symmetrized xx 
electric-field fluctuations by using an obvious analogy of 
relation ^. These electric- field fluctuations are related 
to the XX elements of the Green function tensor via Eq. 
pTj) . The terms proportional to 6{x — x') in Eqs. ((25|) 
and (j27p can be ignored since they originate from the 
short-distance terms proportional to 5{r — r') which do 
not play any role in the large-distance asymptotic. Since 
the combination 



i?x.(0+, 0+) + D,40',0-) - 2D^40+,0-) 
_ AT:h{cqf 1 



OJ Kie2 + K2ei V^i €2 



!a + fi-2 



(28) 



we finally get 



(3h 



PS{uj,q) = ^ coth{(3hLu/2)qHmg{uj,q), (29a) 
An 



'){^,q) 



Ki{uj,q)e2{u}) + K2{uj,q)ei{uj) 
[ei(u;)-62(^)P 



ei(w)e2('^) 



(29b) 



Our task is to determine, for given ei(aj) and £2(^)1 the 
small-q behavior of the function 



f3S{t,q) 



duj 

2^ 



e~"^'f3S{uj,q). 



(30) 



III. ANALYSIS OF ASYMPTOTIC BEHAVIOR 



In any material medium, the complex dielectric func- 
tion e(aj) possesses the symmetry |l9| 



e*{uj) = e{-uj), uj e R. 



(31) 



Denoting e{uj) = e' (uj) + ie" (uj) , where both the real £'{uj) 
and imaginary e"{i-u) parts are real numbers, this means 
that 



where ujp is the plasma frequency defined in Eq. ^ and 
the dissipation constant rj is taken as positive infinitesi- 
mal, 77 ^ 0+ . The Weierstrass theorem reads 

lim — ^ ^V ( -iTi-^Six), X e R, (35) 

r)^0+ x±iri \x J 

where V denotes the Cauchy principal value. It is tempt- 
ing to apply this theorem directly to the representation 
(IM)) . with the result 



e[ico) = l-culr(^-^y e'i[u) 



■KLol-5{uj). (36) 



Although both real and imaginary parts satisfy the nec- 
essary conditions (|32p and ([551) . the expression for the 
imaginary part has no meaning. In the algebraic manip- 
ulations with ei(w), we must therefore keep the positive 
infinitesimal rj in the representation ([34]) up to the end 
and to apply the Weierstrass theorem ([55)) only to the 
final formula. In the vacuum region. 



e2(w) = 1. 



(37) 



To determine the sign of some quantities, we shall some- 
time need the infinitesimal imaginary part of the vacuum 
dielectric constant. To fulfill the required properties ([32]) 
and p3p . the vacuum dielectric constant in fact corre- 
sponds to the limit 



e2{ijj) = 1 + isg-n{LL>)0^ 



(38) 



It is evident that the inverse length k, defined by the 
relation ([24[) . also possesses the symmetry 



K*{uj, q) = k{—oj, q), uj E R. 



(39) 



In terms of the real and imaginary parts, k{uj, q) = 
k'(uj, q) + iK"(uj, q), this symmetry is equivalent to 

K'{u;,q) = K'{-LJ,q), n" {w,q) = -n" {-uj,q). (40) 



e\u) = e'{-u), e"{uj) = -e"{-u:). 



(32) 



The dielectric function has many other general proper- 
ties. Like for instance, for any material medium with 
absorption it holds 



A. Nonretarded limit 



In the nonretarded limit c = cxd, the definition 
becomes ki = K2 = g and Eq. ([29bp reduces to 



e"(w) > 



when cj > 0. 



(33) 



We shall analyze the general fluctuation results of the 
previous section for the model configuration of interest: 
the semi-infinite jcUium in vacuum. As concerns the jel- 
liunr, the dissipation goes to zero in the limit of small 
wave numbe rs pOl . The dielectric function can be thus 
shown to be [21j the Drude one 



ei(w) = 1 



Ijj{lj + irf) ' 



(34) 



5(^,9) 



q 



ei(a;) e2(w) ei(cj) + £2(0;) 



(41) 



Substituting here the dielectric functions ([M)) and (|37p . 
and then applying the Weierstrass theorem ()35p . we ob- 
tain 



Im5(c^,q) = -^sgu{u:)[5{u:^~uj^^)-5{uj^-ul)\ 

= -UJ [2w,,5(w2 _ ^2) _ ^^^(^2 _ ^2)j . (42) 



the frequency of surface plasmons Us is defined in Eq. 
([7]). Using the general formula for the (5- functions 



S[hix)]^J2 



d{x — Xj) 



(43) 



where {xj} are the real roots of h{x), it is a simple task 
to show that relations (|29ap and ([50]) imply the expected 
result (fT2|l . 

The c = oo treatment of our equations is mathe- 
matically accessible and reproduces the previous "non- 
retarded" results. Since the true value of c is finite, 
there exists a limitation for wave numbers q or distances 
A ^ 1/q for which retardation effects are negligible. We 
can find this limitation via a simple dimensional anal- 
ysis of Eqs. (|29p and ([50)1 . The invoked substitution 
f3huj — uj' tells us that there are two dimcnsionlcss quan- 
tities in the theory: 



Phcq, 



f3hujp. 



(44) 



The requirement of the smallness of u ^ 1 is equivalent 
to the natural condition 



Aph <C A, 



(45) 



where Aph oc phc is the thermal de Broglie wavelength 
of photon. The nonretarded limit c —^ oo describes ade- 
quately the region corresponding to the inequality u^ v, 
i.e. 



A< 



OJr, 



(46) 



For rough value of the free electron density in metal 
n ~ lO^^m^^, taking e and m as the charge and mass 



of electron, we have ujp ^ 10 
A < 10"^ 



^^s ^ and so 



lO^d, 



(47) 



where d ^ 10~^m is the mean interparticle distance. We 
see that the distance, over which the retardation effects 
are negligible and so the formula ([12]) is adequate, is rel- 
atively large. 



B. Retarded region 

If we want to deal strictly with the large- distance 
asymptotic behavior of the surface charge correlation 
function, we must take c as a finite number and con- 
sider the region corresponding to the inequality u <ti v, 
or cq -^ Up. This requires a detailed analysis of Eqs. (I^^)) 
and PU)1. 

We start with an algebraic treatment of the function 
g{uj,q) defined by Eq. (|29b[) . Multiplying both numera- 
tor and denominator by Kie2 — K2ei, and using the equal- 
ity 



^l'=2 



^2'^l 



= (£2 -ei) 



<7^(ei +£2) 5-eie2 



(48) 



we obtain 



9 



Kie2 - K2ei 



g2(ei +€2) - (w2/c2)eie2 



(49) 



The denominator of this expression for g is related to the 
surface-plasmon dispersion relation 



q^iei +£2) Teie2 = 0; 



(50) 



for a review, see Ref. f2T!, Considering the metal dielectric 
function ei(w), as defined by Eq. ([M| but for the time 
being with i] = 0, and the vacuum one 62 (w) = 1, Eq. 
(1501) has two solutions 



u;l{q)^icjl/2) + icqr± {ul/2f + [cqf 



1/2 



(51) 



These two dispersion relations are represented, for w > 0, 
in Fig. 2 by solid lines, together with the (dashed) light 
line UJ = cq. The upper solid line w^. always lies above 
the dispersion curve of light in the metal [19| , 



uliq) > ul + {cqf 



(52) 



The lower solid line Lu^{q), corresponding to the disper- 
sion relation of the surface-plasmon polariton, lies always 
below the light line, 



^-(9)<(c<z)' 



(53) 



In the nonretarded limit cq 3> Ws, it approaches the 
nondispersive surface-plasmon frequency uig = ujp/\/2. 
In the retarded region of interest cq <^ uj^, it approaches 
the light line lo = cq. In what follows, we shall need the 
small-g expansions 






{cqf 



[cqY 



0J„ 



{cq? 



{cqf 



(54a) 
(54b) 




FIG. 2: Dispersion (solid) curves u)-{q), uj+{q) and the light 
(dashed) line uj — cq; the plane vector q is represented in units 
of LUp/c, the frequency uj in units of ujp. 



After substituting the metal dielectric function (p4| . 
now completely with 77 being positive infinitesimal, and 
the vacuum dielectric function e2{'-o) = 1 into (j49| . we 
obtain 



where 
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5(w,g) = gi{uj,q) +g2iuj,q), 
KiU}{uj + irf) 



(55) 



(w^ — uPl)[u}'^ — cj^) + iujiyi]' 



92 = {cujpf 



4: 



(oj^ — w^)(w^ — w^) + iujvri 



(56a) 
(56b) 



Here, the functions Ki{Lj,q) and K2(a;,g) are the square 
roots of 



2 _ 2 I ^P ^^ -^P ^V 

'^l ^ 'i ' 2 2 *2 2i2' 

C C C W + 7y 

K2 = g^ - — - I— sgn(w)0+ 



(57a) 
(57b) 



with a positive real part and the new function ^{lo) is 
defined by 

< for w^ = Lo'^, 
v{uj) =u? - 2(cqf { (58) 

> for up' — uj\. 

As indicated, only the sign of this function at the points 
up' = {jp_ and ijp = ijj\ will be needed. 

In order to keep the transparency of algebraic op- 
erations, we introduce the counterparts of quantities 
PS(ijj, q) ([29a|) and ^S{t, q) ^ for each of the compo- 
nents g\ and g2'- 

PSJu},q) = !^coth{(3huj/2)q^lmgj{uj,q), (59) 
47r 



PS,{t,q) 



duj 
2^' 



'f3S,iu;,q), 



j — 1,2. The quantity of interest is 

PS{t,q)^PSi{t,q)+pS2{t,q). 



(60) 



(61) 



The evaluation of the introduced functions is the subject 
of the Appendix. We summarize shortly the obtained 
results in the next paragraph. 

There are two linear in q contributions to f3Si{t,q). 
The first one (jA.5p originates from the discrete level 
of bulk plasmons w = ±Wp, the second one (jA.12p re- 
sults from the integration over the continuous spectrum 
Ci)'^ > ojp + {cqY- It is seen that the two contributions are 
exactly canceled with one another, so that 



As concerns the quantity pS2{t, q), there is only one lin- 
ear in q contribution (|A.20p coming from the integration 
over the continuous spectrum w^ > (cq)"^, i.e. 



f3S2it,q) = ^ 
An 



o{q). 



(63) 



We conclude that the total pS{t, q) has the small-q ex- 
pansion of the classical static type 



m,i) = l^ 



o{q). 



(64) 



In other words, for both static and time-dependent sur- 
face charge correlation functions, the inclusion of retar- 
dation effects causes the quantum prefactor to take its 
universal static classical form. This result holds for any 
temperature. The static t — Q version of (|M|) is clearly 
consistent with the classical finding (O, as it should be 
in the high-temperature limit ph — > ^llj . On the other 
hand, when i 7^ 0, our result ([64)) does not reproduce the 
classical time-dependent one (IT^ , but this feature is not 
against the general principles discussed in the Introduc- 
tion. 



IV. CONCLUSION 

We have studied the long-range decay of the charge 
correlation function on the surface of the conductor in 
vacuum. This problem has been investigated previously, 
for both static and time-dependent correlation functions, 
in the classical and quantum nonretarded regime. Within 
the framework of the fluctuation EM- field theory we have 
shown that the consideration of retardation effects leads, 
for any temperature, to the universal static classical form 
of the asymptotic decay: 



/?^qu(t,|R-R'|) 



1 



1 



|R-R'Hoc 87r2|R-R'|3' 



(65) 



PSi{t,q)^o{q). 



(62) 



independent of t and h. 

As a model system for the conductor, we have used 
the jellium with the simple dispersion relation of Drude 
type ([51)1 . It is not clear at the present stage whether 
the obtained result takes place also for other Coulomb 
or dielectric models with other types of the dispersion 
relation. 

As a byproduct of the formalism, we have obtained a 
very simple formula (|4ip for the g-function, valid in the 
quantum nonretarded regime. This formula reproduces 
the previous microscopic result for the jellium and might 
be used to analyze general Coulomb or dielectric mod- 
els, characterized by their dielectric functions. A general 
analysis might be possible also for the retarded case. 

The extension of the present macroscopic study to 
other domain geometries, e.g. a conductor confined to 
a slab, might be of interest. 

In the near future, we plan to perform a different analy- 
sis of a jellium coupled to the electromagnetic radiation, 
based on the method of collective modes developed in 
Refs. 0,i. 
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APPENDIX 

1. Contributions from gi 

Let us first analyze the contributions of gi{Lj^q) to 
f3Si(t,q), obtained by using Eqs. ((55|) and ([M)) . 

Dividing ni onto its real and imaginary parts, ni = 



k'i + in'l, Eq. ()57ap splits into two relations 



K)^-K')2 = q' + ^^^.^ 



A\^l t^-t 



1(? up' + 77^ 



(A.la) 
(A. lb) 



The function on the rhs of (jA.lbp will help us to choose 
the correct sign of k" which is consistent with the condi- 
tion k!y > 0- We have to distinguish between two cases. 



^1 - "V^P 



{^cqf 



1 < 



u)r] 



k[ 2(? up + r\^ 



(A.2a) 



• u? >ujI + {cqf: 



1 < 



UJJ] 



k'{ 2P Cj2 + 7^2 



>0, 



< = --sgnH^c.2-c.2-(cg)2. (A.2b) 

Here and hereinafter, we adopt the convention that the 
square root of a real positive number has the plus sign 
and neglect all terms of order 77^. Since ujri/{uP + rj'^) oa 
uj6{uj), the integration over uj of the functions which con- 
tain this factor gives zero contribution. 

The function gi, defined by Eq. (j56ap . can be ana- 
lyzed by using the Weicrstrass prescription ((35|) . Since 
Ki(wp) = q and {uj^ - uj^^){u}l ~ ujD = -uj'^{cq)'^, we 
obtain 



Imgi 



—sgn{uj)S{uj'^ - ujp) - [cWpfuj'^K'l 



xV 



up — uj"^ 



V 



V 



w^ — w^ 

(A.3) 

We take into account the lower bound (f52|) for ^^(q) 
and the upper bound (|53p for uP_ (q) and, as before, dis- 
tinguish between two intervals. 



Imgi = -sgn(w) — i^S{u;^-uj;) 



-sgn(w) 



TTCUJpU)'^ 



{ujI-ujI){ujI-u:1) 



X ^ul + (cg)2 - ujH{u:^ - ujI)- (A.4a) 



• U}"^ > ujI + [cqf: 



Ivagi ^ sgn(w) 



CUJpUJ^ 



x^u;^-Lol^{cqrr (^J^^) (A.4b) 

The contributions of the discrete (5-function terms in 
fOa)) to pSi{t,q), defined by Eqs. (PTl and pi)) . are 
easy to find. The bulk-plasmon term in (|A.4a[l . which is 
proportional to 5{uj'^ — Wp), gives a contribution of order 

q- 



in 



cos{ujpt)f{ujp), 



(A.5) 



where the function / is defined by Eq. (|8b|) . The term 
in (|A.4ap . which is proportional to S{uj'^ — uj'^), gives a 
contribution of order q^: 

q^ cWp 

- -— cos(w_i)/(w_)— ^ wr^ T\ 

47r ' {lj± - uj^){ujl ~ ujt) 



{cqf 
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q c 

47r UJr, 



+ o(g2). (A.6) 



To find the contribution of the term in (jA.4bp . which is 
continuous in w, is a more complicated task. In the small- 

q limit, two functions l/(w^ — ujT) and \ ui"^ — uj'i — (cg)^ 

are close to their singular points just at the lower bor- 
der of the integration uj'^ — uj'^ + (cq)'^. Moreover, the 
principal value 7^(1/(^2 — w^)) has also to be taken with 
caution; see the small-g expansion (j54bp of w^. Based 
on the above information, we split the whole interval of 
w-values uj'^ > ujp + {cq)^ onto the small one symmetric 
with respect to w^, 



h : ujj^—l<uj <uj'^+l, I = uj'^ — [ujp + {cq) 



and the infinite one 



uj^+l< up. 



UJp 

(A.7) 



(A., 



We first perform the integration of the three "problem- 
atic" quickly changing functions over the interval /i . In- 
troducing the variable u via w^ = uj'i^-\-u, the contribution 



from the w-integration over the interval /i is the factor 
(cq"^ /A'k'^) cos{ujpt){(3hujp/2) coth{(3hu;p/2) multiphed by 



The function 52 is defined by Eq. (|56bp . Using the 
Weierstrass theorem, we obtain 



du 7 rT5 

I u + {cqY + I 



Vu + l- 



,2 _|_ ^2 



{cq? 
[cqY 



du 



u \ 



vr 



\/2-sinh"^(i; 



o(i). 



(A.9) 



This means that the contribution of (|A.4bp to j3Si{t,q) 
coming from the interval /i is of order q^ . The integration 
of the term (jA.4bp over the interval I2 can be represented, 
with the substitution uj"^ = lo"^ + 1 + v and in the small-g 
limit, as follows 






dv cos 



[t^pl^v 



y/v + 2l 



ijji+v 



X coth I —^Juj^ + v 



V + I V + {cqY 



(A.IO) 



Performing the next substitution v — [cq^w, considering 
the g — > limit and evaluating the integral 



dw—= 

V"^ w + 1 



(A.ll) 



the total contribution of the term (jA.4b[) to f3Si (i, q) 
reads 



— COs{ujpt)f{ujp) +o{q). 



(A.12) 



The linear in q part of this contribution is exactly can- 
celed by the bulk-plasmon contribution (|A.5|) . We there- 
fore conclude that 



PSiit,q)^o{q). 



(A.13) 



lmg2 = {aVpY^V { . ^ A v 



,2- ,,,2 



2 _„,2 



-TT{cUJp)'^K2Sgn{uJl')S[{Lj'^ — Uj'^){uj'^ — ^+)]- 

(A.15) 



We thus have 



• uj"^ < (cq)^: 



(A.16a) 



• uj"^ > (cq)^: 



lmg2 = -cujI^/u!^ - {cqysgn{uj) 



xV 



9 2 / I 9 2 



.(A. 16b) 



The contribution of the term (|A.16a|) to PS2 {t, q) reads 

-cos(c.-0/(a;-) ^2(^^„^2) 
t 

All (jJr, 



X^{cqr-Lul = ^- + Oiq'). (A.17) 



The contribution of the term (|A.16b[) can be represented, 
after the substitution u = uj'^ and in the small-g limit, as 
follows 



phq^ 



i-ccol) 



StT^ P J, y- Vu 



du 



3(tVw) coth(/3ftV^/2) 



xV \ ^—] ^^==.{k.l8) 
u-uj^J ^u - {cqY 



2. Contributions from 32 



Making the next substitution u — [cqY'v and evaluating 
the integral 



We proceed by a relatively simple analysis of the con- 
tributions of .g2((j-', q) to /3S2(t, q). 

Writing ^2 = ^2-!- iK2 , we have to distinguish between 
two cases. 



• w^ < (cq)^: 



4 = -v/R^^^, 4-0; 



• w^ > (cq)^: 



(A.14a) 



4^0+, 4 = ~-sgn{uj)y^uj^ - (cqy. (A.14b) 



dv- = TV. (A.19) 

1 V ^Jv -I 



the contribution of the term (|A.16bp to (3S2 {t, q) is found 
in the classical static form 



47r 



oiq). 



(A.20) 



In view of Eqs. (|A.17p and (|A.20|) we conclude that 



(3S2it,q)^^ + oiq). 



(A.21) 
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